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Abstract 

The decays B° — > D^jK~ and B° — > Djjir + are studied for the first time. A significant sig- 
nal is observed in the B° -» Z?; j (2317) + K- decay channel with 5(5° -» D* J (2317) + K _ ) x 
H( J D* J (2317) + -» D+vr ) = (5.3±i;| ± 0.7 ± 1.4) x 1(T 5 . No signals are observed in the 
W -» D* i7 (2317)-7r + , B° -» D sJ (2460)+K-, and 5° -► D sJ (2460)-tt + decay modes, and up- 
per limits are obtained. The analysis is based on a data set of 140 fb _1 collected by the Belle 
experiment at the asymmetric e + e~ collider KEKB. 

PACS numbers: 13.25.Hw, 14.40.Nd 
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Two narrow resonances denoted as D!r(2317V t " and _D s j(2460) + have been observed re- 
cently in e + e~ continuum interactions 0, S 0| • These resonances were initially seen in 
the D*j(2317)+ -> D+vr , D sJ (2460)+ -> D+7 and L> sJ (2460)+ -> D* + 7r° decay modes 
0, and their quantum numbers were tentatively classified as J p = + for D*j(2317) + and 
J p = 1 + for D s j(2AQ0) + . However, the measured masses are significantly lower than the 
values predicted within potential models for + and 1 + states p. The D s j mesons were 
also observed in B — > DD s j decay modes with branching fractions an order of magnitude 
less than those for B — ► DD S decay modes with a pseudoscalar D s 0]. Angular analysis 
of B — > DD s j(2A60) + favors a spin 1 assignment for D s j(2A60) + . There has been a signif- 
icant effort to explain the surpising D s j masses J6|L and some authors have discussed the 
possibility of four-quark content in the Dfj |H 0, 0, [ill E3 ■ 

In this Letter we report the results from a search for B° — > DfjK~ and B° — > Djjir + 
decays, where D^j mesons are reconstructed in the modes -D*j(2317) + — > D+n and 
D s j(2A60) + — > Dfj. Measurements of the corresponding decays B° — > D+K~ and 
B° — > D~ s ~n + have been reported recently by Belle 0] and BaBar [l4^ . 

The decay mode B° — > D^ JS) 7i + can be described by a 6 — > w tree diagram. Within the 

factorization approach the branching fraction ratio R w +/o+ = B(B° — > D~n + )/B(B° — > 
D~D + ) is predicted to be (0.424 ± 0.041) • (V^/V^l 2 and can be used to obtain the ratio 
of Cabbibo-Kobayashi-Maskawa matrix elements Assuming similar ratios R k +/d+ 

for D s and D s j mesons, only a few B° — > Djjir + events would be observed in the current 
Belle data sample. 

The decays B° — > D~^,j^K~ are of special interest because the quark content of the initial 
B° meson (bd) is completely different from that of the D^^K~ final state (cssu), indicating 
an unusual configuration with both initial quarks involved in the weak decay. Branching 
fractions for the pseudoscalar meson B(B° — » Dt K~) = (4.6+j' 2 ± 1.3) • 10~ 5 and 
(3.2 ± 1.0 ± 1.0) • 10~ 5 were measured by the Belle jl3j and BaBar [14J collaborations, re- 
spectively. Predictions for this branching fraction have been obtained assu ming a dominant 
contribution from a PQCD factorization W exch ang e process (Fig. QJi) 0, Il7f or, alterna- 
tively, from final state interactions (Fig. [TJd) jl8L llflj . and cover the range from a few units 
times 10~ 6 to 10~ 4 . If the D s j mesons have a four-quark component then the tree diagram 
with ss pair creation (shown in Fig. [TJ:) may also contribute. 




FIG. 1: Diagrams describing B — * D^jK decay. 

The analysis was performed using a 140fb _1 data sample containing (152.0 ± 0.7) x 10 6 
BB pairs. The data were collected with the Belle detector at KEKB |20(, an asymmetric 
energy double storage ring collider with 8 GeV electrons and 3.5 GeV positrons. Belle 
is a general-purpose large-solid-angle detector that consists of a three-layer Silicon Vertex 
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Detector (SVD), a 50-layer Central Drift Chamber (CDC), an array of Aerogel Cerenkov 
Counters (ACC), a Time of Flight counter system (TOF), and a CsI(Tl) Electromagnetic 
Calorimeter (ECL) located inside a superconducting solenoidal coil with a 1.5 T magnetic 
field. The detector is described in detail elsewhere \2l\ . 

Charged tracks are required to have momentum p > 100 MeV/c 22l | and impact pa- 
rameters less than 2 cm radially and 5 cm in the z direction |23| with respect to the 
interaction point. Kaon and pion mass hypotheses are assigned using a likelihood ratio 
£>k/tt = £k/(£k + £tt), obtained by combining information from the CDC (dE/dx), ACC, 
and TOF systems. We require Cr/tt > 0.6 (Ck/k < 0.6) for kaon (pion) candidates jUf. 
With these requirements the identification efficiency for particles used in this analysis varies 
from 91% to 86% for kaons and from 98% to 94% for pions, decreasing as the momentum 
increases. 

ECL clusters with a photonlike shape and energy larger than 50 MeV, that are not 
associated with charged tracks, are accepted as photon candidates. Photon pairs of invariant 
mass within ±12 MeV/c 2 (~ 3a in the ir° mass resolution) of the tt° mass are considered 7r° 
candidates; the 7r° momentum is required to be larger than lOOMeV/c. 

Kg candidates are formed from tt + it~ pairs with an invariant mass within ±10MeV/c 2 
(~3o~) of the nominal Kg mass and a common vertex displaced from the interaction point 
by more than 0.2 cm in the plane perpendicular to the beam direction. A common vertex for 
the two tracks in the plane perpendicular to the beam direction was found; the difference 



in z coordinates of the measured pion tracks at this point was required to be less than 



2 cm. The angle a between the Kg flight and momentum directions is required to satisfy 



cos a > 0.8. 

Invariant masses of K*° —>■ K + tt~ candidates are required to be within ±50MeV/c 2 of 
the nominal K*° mass; those of <ft ~ ► K + K~ candidates, within ±12MeV/c 2 of the <fi mass. 
.D+ mesons are reconstructed in the 07r + , K*°K + and K$K + decay channels; a mass window 
of ±12MeV/c 2 (~ 2.5cr) is imposed in each case. 

The D s j mesons are reconstructed in the _D*j(2317) + — > D+n° and _D s j(2460) + — > D+j 
decay modes. To select a -D*j(2317) + , the candidate mass difference AM (Z}*j(2317) + ) = 
M(L>+7T°) - M(D+) is required to lie within ±20MeV/c 2 of 348.6 MeV/c 2 (~ 3.0a). To 
select a £> sJ (2460)+, we require AM(D sJ (2460)+) = M(D+-f) - M(D+) within ±30MeV/c 2 
of 487.9 MeV/c 2 (~ 2.5cr). The D s j and D s mass differences were taken from [^]. 

We then form B° — > DfjK~ and D~jH + candidates and extract the signal using the energy 
difference AE = E% M - Eg^ m and beam-constrained mass M bc = ^/(Egf m ) 2 - (pg M ) 2 ; 
E% M and p c B M are the energy and momentum of the B candidate in the center-of-mass (CM) 
system and Eg^ m is the CM beam energy. Only events within the intervals M^ c > 5.2 GeV/ c 2 
and \AE\ < 0.2 GeV are used in this analysis. The B meson signal region is defined by 
\AE\ < 0.04 GeV and 5.272 GeV/c 2 < M bc < 5.288 GeV/c 2 . 

Combinatorial background for channels involving the _D s j (2460) + was further suppressed 
by requiring cos#d s7 < 0.7. The helicity angle 6*d s7 is defined as the angle between the di- 
rection opposite the B momentum and the Df momentum in the rest frame. This re- 
quirement rejects 49% of background events and only 6% of signal events, assuming J p = 1 + 
for the _D s j(2460) + . The uncertainty due to this assumption is included in the systematic 
error. 

For events with two or more B candidates, the Df and 7r° candidates with invariant 
masses closest to their nominal values and the B daughter K + or n~ candidate with the 
best Ck/tt value are chosen. With these requirements no multiple entries are allowed for the 
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FIG. 2: AE (left) and AM(D sJ ) (right) distributions for 5° decays to (a) D*j(2317) + K~ , 
(b) 0*^(2317)-^+, (c) D sJ (2i60) + K- and (d) £> sJ (2460)-7r+. Tight requirements on M bc , and 
AM(D s j) (left) or AE 1 (right) are applied; see the text. 



_D*j(2317) + channels and, according to the Monte Carlo (MC) simulation, less than 1% of 
selected events will have two B candidates in channels with D s j(2A60) + . No multiple entries 
are found in the data. 

After this selection the principal background is from e + e~ — > qq continuum events 
(q = u,d,s, or c). We exploit the event topology to separate BB events (spherical) from 
the continuum background (jetlike). The ratio of the second and zeroth Fox- Wolfram mo- 
ments |2jJ of all particles in the event is required to be less than 0.5. For such events, we 
form a Fisher discriminant from six modified Fox- Wolfram moments jHI]. A signal (back- 
ground) likelihood Cs (£bg) is obtained using signal MC (sideband) data from the product of 
probability density functions for the Fisher discriminant and cos 6 b, where 9 b is the B flight 
direction in the CM system with respect to the z axis. We require 1Z = Cs/ (Cs + Cbg) > 0.4 
for Df — > K*°K + and TZ > 0.25 for the other decay modes, which have lower back- 
grounds. For the B° — > D*j(2317) + K~ mode these requirements retain 92%, 85%, and 95% 
of signal events while removing 47%, 67%, and 64% of continuum events, for — > 07r + , 
K*°K + , and K S K+, respectively. The fract ions retained (or removed) for the other B decay 
modes are similar, varying by a few percent. 

The AE and AM(D s j) distributions for the various DfjK~ and D~jir + combinations are 
shown in Fig.Elfor the range 5.272 GeV/c 2 < M bc < 5.288 GeV/c 2 . To obtain the AM(D sJ ) 
distributions we relax the AM(D s j) requirements and apply a tight selection on AE. Each 
AE distribution is fitted by a Gaussian with zero mean and a width fixed from MC data 
to describe the signal, and a linear background function. The AM(D s j) distributions are 
described by signal Gaussians with widths fixed from MC data and mass differences fixed 
to 348.6 MeV/c 2 or 487.9 MeV/c 2 , and linear backgrounds. A clear B° D* sJ (2317) + K~ 
signal is observed; no significant signals are observed in the remaining modes (Fig. |2J). The 
B° yields, based on fits to histograms combining all three Df decay modes, are listed in the 
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TABLE I: Signal yields, efficiencies, product branching fractions (or limits), and significances for 
the B° — > D^jK~ and D~j-k + decay modes. Only statistical errors are shown. Product branching 
fractions are obtained from AM{D s j) fits: see the text. 
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W -» D* t/ (2317)-7r+ 


9 q+3.3 


0.5 ±3.3 


27.6 ±1.3 


< 2.5(90%C.L.) 




B° -» D sJ (2460)+i^- 


2-oiil 


1.0 ±2.9 


56.5 ±2.4 


< 0.94 (90%C.L.) 




S° -> D sJ (2460)"7T + 


-i-9±S:S 


-3.9 ±4.1 


65.6 ±2.6 


< 0.40(90%C.L.) 





last four lines of Table [H 

Various studies are performed to confirm the B° — > D*j(23l7) + K~ signal. The parame- 
ters of the signal peak are allowed to float in the AM(D s j) fit: a mean (351.2 ± 1.6) MeV/ c 2 
and a width (6.0 ± 1.2) MeV/c 2 are obtained, in good agreement with the MC expectations, 
(348.5 ± 0.3)MeV/c 2 and (6.1 ± 0.2)MeV/c 2 . (The mass of 2317.5 MeV/c 2 and the zero 
width are used in the MC simulation of the _D*j(2317) + signal.) Good agreement is also 
obtained for the signal position and width in AE and M^ c . 

To check for a possible background contribution due to a random combination 
of a -D*j(2317) + meson and a kaon, AM(D*j(2317) + ) distributions are obtained for 
events from the AE sideband 0.05 GeV < \ AE\ < 0.2 GeV, and the M bc sideband 
5.2GeV/c 2 < Mb c < 5.26GeV/c 2 . After rescaling the fit results to the B signal region, 
the background contribution is estimated to be —0.8 ± 0.7 (1.1 ± 0.9) events using the AE 
(Mb c ) sideband. As both AE and M bc requirements are applied in the AM(D*j(2317)) fit, 
this background contribution is estimated to be less than one event, and treated as a source 
of systematic error. 

The shape of the background-subtracted and efficiency-corrected cos#D sVr distribution for 
B° — ► D*j(2317) + K~ decay is compared with those predicted for possible _D*j(2317) + quan- 
tum number hypotheses. (The helicity angle 9r> g7T is denned as for #d s7 , with n° substituted 
for 7.) The distribution is expected to be flat if the _D*j(2317) + has J p = + , or to have the 
form cos 2 9d s tt m the 1~ case; within large errors, it is consistent with a constant. A fit gives 
X 2 = 1-44 for a constant and \ 2 = 4.72 for cos 2 8d s7T , for four degrees of freedom. A larger 
data sample is required for a statistically significant separation of the two hypotheses. 

Signal yields, efficiencies, branching fractions and significances for the studied decay chan- 
nels are shown in Table HI The B° D*j(23l7) + K~ branching fraction is obtained using 
a simultaneous fit to the AM(D*j(2317)) distributions for the three decay channels, 
with independent background descriptions, but common values for the signal width (fixed 
from MC) and peak position (allowed to float). The branching fraction thus obtained is in 
good agreement with the values from the AE and Mb c fits. Efficiencies include all interme- 
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diate resonance branching fractions |26j and were obtained from MC simulation, assuming 
J p = 0+ for the D*j(2317) and J p = 1 + for the D sJ (2A60). We assume equal production of 
neutral and charged B mesons. The significance is defined as a/— 2 ln(£ /£max), where £ max 
and Cq are likelihoods (corrected for the number of degrees of freedom) for the best fit and 
zero signal yields, respectively. The upper limits are obtained using fits to AM(D s j) distri- 
butions, with fixed signal positions and widths. We use the Feldman- Cousins method |2^ |. 
assuming a Gaussian distribution for the statistical error. The upper limit is then increased 
by 29% (the sum in quadrature of the experimental systematic error and the uncertainty in 
the .D+ branching fraction scale). The systematic error is treated in a conservative way in 
order to avoid Bayesian assumptions about its probability distributions. Other methods for 
upper limit determinations agree with the values obtained here within 5% for the first two 
upper limits (that have positive signals) and within 15% for the last upper limit. 

The main result of this study is the measurement of the product branching fraction 
B(B° -> D* sJ (2317) + K-) x B(D* aJ (2317)+ -> D+tt°) = (5.3^| ± 0.7 ± 1.4) • 10~ 5 . The three 
error terms are the statistical uncertainty, the total systematic error, and the uncertainty 
due to branching fractions; this last term is dominated by the ~ 25% uncertainty in 

The major sources contributing to the systematic error are uncertainties in efficien- 
cies of charged track reconstruction (1% x iV trac ks), particle identification for charged pi- 
ons (2% x N n ) and kaons (2-3% x Nk)', the photon and 7r° reconstruction efficiencies and 
energy scale (5%); the K s vertex reconstruction (3%); the efficiency of the AE (2%) and 
topological likelihood ratio (TZ) selections (3%); the background and signal shape definition 
for the B signal (3%); the background subtraction (6%); the change in reconstruction effi- 
ciency for the different Dfj quantum number assumptions (4%); the statistical uncertainty 
of the MC sample used to determine efficiency (4%); and the uncertainty on the number of 
BB pairs (0.5%). These uncertainties were added in quadrature to obtain a total systematic 
error of 13%. 

In conclusion, the B° — > D^jK~ and D~j7t + decay modes were studied for the first 
time. The B° — > D*j(2317) + K~ mode was observed, with a product branching frac- 
tion B(B° -> D;j(2317)+K-) x B(D* sJ (23l7) + -> D+n°) = (5.3^ j ± 0.7 ± 1.4) • 10" 5 . Re- 
cent measurements imply that the _D*j(2317) + — > D+tt° channel is dominant and the 
_D s j(2460) + — > fraction is around 30%. Taking into account these approximate val- 
ues, we can conclude that B(B° — > D* sJ {23\7) + K~) is of the same order of magnitude as 
B(B° — > D+K~) and at least a factor of two larger than the B° — > D s j(24Q0) + K~ branch- 
ing fraction, in contrast to the naive expectation that decays with the same spin-doublet 
D*j (231 7) + and _D s j (2460) + mesons would have similar rates. No significant signals for 
B° — > D~j7f + decays were seen. 
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